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Abstract 

We report the first observation of the radiative transition hb(lP) — > rjb(lS)j, where the hb(lP) 
is produced in Y(55) — > hb(lP)ir + Tr~ dipion transitions. We measure the 7#>(15) mass to be 
(9401.0 ± 1.9±2:t) MeV/c 2 with a width of (12.4j;^t^ 4 5 ) MeV and a decay branching fraction of 
B[hb(lP) ->■ rj b (lS)-f] = (49.8 ± 6.8l5° 2 9 )%. The measured 776(15) mass corresponds to a hyperfine 
splitting of (59.3±1.9lf 4) MeV/c 2 . This value deviates significantly from the current world average 
obtained from measurements of T(35) — > rjb(lS)j and T(25) — > 775(15)7 reactions. We also 
report updated results for the hb(lP) mass (9899.0 ± 0.4 ± 1.0) MeV/c 2 and its hyperfine splitting 
(0.8 ± l.l)MeV/c 2 . These measurements are performed using a 121.4fb _1 data sample collected 
at the peak of the T(55) resonance with the Belle detector at the KEKB asymmetric-energy e + e~ 
collider. 

PACS numbers: 14.40.Pq, 13.25. Gv, 12.39.Pn 
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INTRODUCTION 



Recently Belle reported the first observation of the /i&(1P) and hb(2P) states [11]. The 
radiative transition to the r] b (lS) is expected to be one of the dominant decay modes of 
the h b (lP): ggg/r/^lS)^ /gg^f = 57/41/2% [2]. Belle's large h b (lP) sample provide an 
opportunity to study the rjb(lS), which is the ground state of the bottomonium system with 
bb spin and orbital angular momentum equal to zero. The hyperfine splitting defined as 
AMhf[?7&(1'S')] = M[T(15')] — M[rjb(lS)\ provides a test of spin-spin interactions [3J. The 



rjbilS) was first observed by BaBar |4j,|5] and confirmed by CLEO jfj]. Its mass is found to 
be higher than theoretical predictions [7, The tension between experimental results and 
predictions strongly motivates further experimental studies of the rjb{lS). We note that no 
experimental information is available on the rjb(lS) width. 

We report the first observation of the radiative transition hb(lP) — > 7^(1 5)7 and mea- 
surements of the r}b(lS) mass, width and decay branching fraction. We use a 121.4 fb _1 
data sample collected at the peak of the T(5S) resonance (\/s ~ 10.865 GeV) with the Belle 
detector [9( at the KEKB asymmetric-energy e + e~ collider 10]. 



METHOD 

In the decay chain T(5S) —> Z^n^ , — > hb(lP)ir + , hb(lP) —> r}b(lS)j we reconstruct 
only the tt~, 7r + and 7. Here Z^ denotes the Zft(10610) + and _Zj,(10650) + , the two charged 
bottomonium-like resonances first identified by Belle in Ref. In this decay the typical 
momenta of the 7r _ , tt + and 7 are 240MeV/c, 730MeV/c and 500MeV/c, respectively. 

We define the missing mass of X (where X = 7r + 7r~ or 7r + 7r~7) as 

M miss (X) = y/(E c . m .-E* x ) 2 -p%, (1) 

where E c m is the center-of-mass (cm.) energy and E* x and p* x are the energy and momen- 
tum of the X system measured in the cm. frame. 

The 7T + 7r _ 7 combinations from the signal decay chain form a cluster in the M m i SS (7r + 7r~7) 
versus M miss (7r + 7r~) plane centered at M[rjb(lS)] and M[hb(lP)], respectively (see Fig. CQ. 
In this plane there is a vertical band due to correctly reconstructed 7r + 7r~ combinations 
and misreconstructed 7's, and a slanted band due to correctly reconstructed 7's and mis- 
reconstructed 7r + 7r~. We introduce a new variable AM miss (7r + 7r~7) = M miss (7r + 7T~7) — 
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M m i SS (7r + 7r~) + m[hb(lP)]. The advantage of the new variable is that the band with cor- 
rectly reconstructed 7's and misreconstructed 7t + tt~ combinations becomes horizontal (see 
Fig. [Q and the variables AM miss (7r + 7r~7) and M miss (7r + 7r~) are not correlated for signal 
events. 





M mis tn + n), GeV/c' 



9.8 9.85 9.9 9.95 10 



M miss (7rV), GeV/c 2 



FIG. 1: Results of Monte- Carlo simulation for signal hb(lP) — > ^(IS)"/ transitions. M m i ss (7r + 7r~7) 
vs. M m i ss (7r + 7r~) distribution (left) and AM m i ss (7r + 7r~7) vs. -/Vf m i ss (7r + 7r~) distribution (right) for 
all 7r + 7r~7 combinations in the event. 



It may be possible to perform a two dimensional fit to the AM miss (7r + 7r~7) vs. 
M m i SS ('K + 'n-~) distribution. However, we follow a more intuitive approach. We divide the 
AM miss (7r + 7r~7) vs. M miss (7r + 7r~) plane into 10MeV/c 2 wide horizontal slices, project each 



slice onto the M„ 



77 r 7T 



axis and fit the M„ 



7T 1 7T 



distribution. We thus find the 



dependence of the hb(lP) yield on the AM m i ss (7r + 7r 7) variable. We then search for the 
rjb(lS) signal as a peak in this distribution. 



SELECTION 



We start with hadron event selection 



12| . Continuum e + e — > qq {q = u, d, s, c) 



background is suppressed by requiring that the ratio of the second to zeroth Fox- Wolfram 
moments satisfies R2 < 0.3 [3]. The 7r + 7r _ selection requirements are the same as in the 
hb(lP) and hb(2P) observation paper We consider all positively identified 7r + 7r~ pairs 
that originate from the vicinity of the interaction point. We require that the hh{^-P) is 



produced via an intermediate Z b , 10.59 MeV/c 2 < M miss (7r) < 10.67 MeV/c 2 This 
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requirement significantly reduces the background (by a factor of 5.2) without any significant 
loss of the signal. We consider all 7 candidates, and apply a ir° veto, |M(77 2 ) — m^o\ > 
13MeV/c 2 with E r2 > 75 MeV, where 72 is any photon candidate in the event. The values 
of the cuts were optimized based on the figure of merit FoM = -f=, where S is the number 

V B 

of signal events in the signal Monte-Carlo (MC), B is the number of background events 
estimated from a small fraction (0.1%) of data. 

STUDY OF INCLUSIVE h b (lP) SIGNAL 

The fit to the inclusive M m i SS (7r + 7r^) spectrum (before combining tt + it~ and 7 candidates) 
with the requirement of the intermediate Z b is shown in Fig. |2j We use the same fit procedure 



x 10 • 




10.1 



M miss (7tV), GeV/c' 



FIG. 2: M m i ss (7r + 7r ) spectrum with the requirement of an intermediate Zf, (top) and residuals 
(bottom). 

as in Ref. [lj]. The fit function consists of four components: the h b (lP) signal, the T(2S) 



S 



signal, a reflection from the T(3S) — > T(1S)tt + tt~ decay, and combinatorial background. 
The shapes of the peaking components are determined from the analysis of the /i + yU - 7r + 7r~ 
data sample, that contains the T(nS) — > (n = 1, 2, 3) decays. The signals are found to 

have tails that account for about 8% of the yield and are due to the initial state radiation of 
soft photons. The hb(lP) and T(2S) intrinsic widths are negligible compared to the detector 
resolution, therefore the signals are described by a Crystal Ball function with width a = 
6.5MeV/c 2 and 6.8MeV/c 2 , respectively. The width (cr) of the hb(lP) signal is determined 
from linear interpolation in mass from the widths of the T(nS) peaks. The tail parameters 
of the hb(lP) signal are assumed to be the same as for the T(2S) signal. The T(3S) — > 
Y(lS') 7r+7r_ reflection is described by a single Gaussian function with a width of a = 18 MeV. 
The combinatorial background is parameterized by a third order Chebyshev polynomial. We 
perform binned \ 2 At using 1 MeV/c 2 bins, though for clarity we display the data in 5 MeV/c 2 
bins. The results of the fit are given in Table [H The confidence level of this fit is 56%. 

TABLE I: The yield and mass of the hb(lP) from the fit to the inclusive M m i ss (7r + 7r~) spectrum. 





Yield, 10 3 


Mass, MeV/c 2 


h b (lp) 

T(3S) T(15) 
T(2S) 


(61.3 ±3.1^|) x 10 3 
(13 ± 7) x 10 3 
(54.8 ± 3.9) x 10 3 


(9899.0 ±0.4 ±1.0) MeV/c 2 
9973.0 MeV/c 2 
(10021.1 ±0.5) MeV/c 2 



To estimate systematic uncertainty on the hb(lP) parameters we vary the Chebyshev 
polynomial order (±1, ±2); and fit range (we reduce it to 9.98MeV/c 2 and exclude all 
peaking components except for the hb(lP) signal). We also introduce a correction factor for 
the signal width and allow it to float (we find / = 0.99 ± 0.07). We use a signal tail shape 
not only from the T(2S) (the default case), but also from the T(IS') and T(3S). A summary 
of the systematic uncertainties is given in Table [TTJ For the mass measurement we introduce 
an additional ±1 MeV/c 2 uncertainty due to possible local variations of background shape 
as estimated in Ref. |lj using deviations of reference channels from the PDG values. The 
new value for the hb(lP) mass corresponds to a hyperfine splitting AMn F [h b (lP)] = (0.8 ± 
1.1) MeV/c 2 , where statistical and systematic uncertainties are added in quadrature. 
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TABLE II: Systematic uncertainties in the h b (lP) parameters from various sources. 







Polynomial 


Fit 


Signal 






order 


range 


shape 


N[hb(lP)], 


10 3 


+ 1.8 
-0 


+1.1 

-0 


+0.5 
-0.3 


M[h b (lP)} 


MeV/c 2 


+0.1 
-0 


+0.1 
-0 


+0 
-0.1 



EXTRACTION OF rj b (lS) SIGNAL 

To extract the rj b (lS) signal we fit the M miss (7r + 7r~) spectra in the AM miss (7r + 7r~7) bins. 
In the fit function we fix the masses of signals at the values given in Table HI We use a 
10MeV/c 2 bin size in AM miss (71+71-7). The results for the h b (lP), T(2S) and T(35) ->■ 
T(15') reflection yields as a function of the AM miss (71+71-7) are shown in Fig. [3j The h b {lP) 
distribution shows a clear peak at 9.4GeV/c 2 that we identify as the r] b (lS) signal, while 
the other distributions do not exhibit significant structures. 

We search for peaking backgrounds in the AM miss (71+71^7) distribution of the h b (lP) 
yield. We use a MC simulation for generic T(55) decays and consider separately h b (lP) — > 
ggg, h b {lP) —> ggj and e+e~ —> jisrT(3S) processes. We do not find any sources of peaking 
background. 

CALIBRATION 

We use the B + — > xaK + , Xci — > J/ipJ, J/ip — > e+e~ / fi + /i - data sample for calibration. 
We require that the kaon and lepton candidates are positively identified and originate from 
the vicinity of the IP. For the J/ip — > e + e~ mode we attempt to reconstruct and recover 
bremstrahlung photons. The mass window around the nominal J /if) mass is ±30MeV/c 2 
(±50 MeV/c 2 ) for the fi + n~ (e + e~) mode. We perform a mass constrained fit to the J/ip 
and Xci candidates. We require \AE\ < 30 MeV and M bc > 5.27MeV/c 2 , which are loose 
requirements. The AE sidebands are defined as 40 < \AE\ < 100 MeV. The background is 
efficiently suppressed by the requirement cos 6 1 > —0.2, where 9 1 is the helicity angle of the 
Xd defined as the angle between the 7 momentum and Xci boost direction in the Xci rest 
frame. The particular value of the cut is chosen so that the average 7 energy in the cm. 
frame is 500 MeV/c 2 , i.e. is equal to the average energy of the photon in the h b (lP) — > r) b (lS) 
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AM miss( 7l+7I », GeV/c 2 

FIG. 3: The results for the h b {lP) (top), T{2S) (middle) and T(35) -> T(15) reflection (bottom) 
yields as a function of AM m i ss (7r + 7r _ 7). 

transition. The AE distribution with the requirement 3.44 < M(J/i/)j) < 3.54GeV/c 2 and 
M(J /ipl) distribution for the AE signal region and for the normalized sidebands (see Fig. H]) 
indicate, that the background is small. 

We parameterize the J / ^7 mass distribution in MC by a Crystal Ball function [ijj with an 
asymmetric core and with tails on both sides (8 parameters in total). This parameterization 
describes MC reasonably well (see Fig. [5]). For the core we find o\ = 11.9MeV/c 2 (left side) 
and 02 = 6.8MeV/c 2 (right side). 

We fit the M(J /ipj) distribution in data fixing the o\ and 02 parameters and introducing 
a shift of the peak position and width correction factor (see Fig. We find for the shift: 
-0.7 ± 0.3±o;4 MeV/c 2 and for the width correction factor: 1.15 ± 0.06lo:o6- The systematic 
uncertainty is estimated (1) by varying the fit interval; (2) by using only the left normalized 
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FIG. 4: (Left) AE distribution for the selected B + — > XciK + decay candidates; signal and 
sidebands regions are hatched. (Right) M(Jfip'y) distribution for the AE signal region (points 
with error bars) and for the normalized sidebands (open blue histogram). 



8000 



6000 



4000 



2000 - 




44 3.46 3.48 3.5 3.52 3.54 
M(J/ V y), GeV/c 2 



FIG. 5: The M(J/ipj) distribution in MC (left) and data AE signal region with AE sidebands 
subtracted (right). The fit is described in the text. 



AE sideband for subtraction or only the right sideband; (3) by varying the oyjo^ ratio in 
the parameterization: we use the 01/02 ratio from the fit to the distribution in (i) data and 
(ii) in the h b (lP) %(15) 7 MC. 
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776(15) MASS AND WIDTH 



We fit the AM miss (7r + 7r 7) distribution to a sum of an rjb(lS) signal component and a 
combinatorial background contribution (see Fig. E]). The signal is a non-relativistic Breit- 




AM miss (7cVY), GeV/c 2 

FIG. 6: AM m i ss (-7r + 7r~7) distribution of the hb(lP) yield with fit result superimposed. 

Wigner (this parameterization is chosen to simplify comparison with BaBar and CLEO 
results) convolved with the calibrated resolution function (we use the shift and width cor- 
rection factor determined from the xa signal in data). The combinatorial background is 
parameterized by an exponential function. We verify the fit procedure using MC. The fit 
results are shown in Table llHl The confidence level of the fit is 77%. The significance of the 

TABLE III: The yield, mass and width of the rjb(lS) signal from the fit to the AM m i ss (-7r + 7r~7) 
distribution of the hb(lP) yield. 



N[v b (lS)] 


(21.9 ± 2.0±f •?) x 10 3 


M[ Vb (lS)} 


(9401.0 ± 1.9±H) MeV/c 2 


r[vb(is)} 


(12.4l^+^ 5 ) MeV 



T]b(lS) signal is 14 a. 

To estimate the systematic uncertainty we vary the fit range in the fits to the M m i SS (7r + 7r~) 
distributions (instead of the default range 9.8 — 10.1 GeV/c 2 we use 9.8 — 9.98 GeV/c 2 ); we 
also vary the polynomial order in these fits (we increase the order by one and by two); we 
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vary the binning of the AM m i SS (7r + 7r~7) distribution (we scan the starting point of bin with 
1 MeV/c 2 steps); we also vary the fit range in the fits to the AM m i SS (7r + 7r~7) distribution; we 
vary the parameterization of the combinatorial background in the fits to the AM miss (7r + 7r~7) 
distribution (in addition to a single exponential we use the sum of two exponentials, second 
and third order polynomials, we also use the function exp(^^ =0 CkX k ), i.e. an exponential 
of a polynomial, with n = 2,3,4); to take into account the uncertainty in the resolution 
function we vary the width correction factor by ±1 o (we combine its statistical and system- 
atic uncertainty in quadrature) and we take into account the uncertainty in the mass shift 
parameter; we vary the selection criteria [R 2 cuts: 0.25, 0.3 (default), 0.35, 0.4; 7r° veto, cut 
on |M( 7 7 2 ) - M w o|: 10, 13 (default), 15MeV/c 2 , cut on E j2 : 50, 75 (default), 100 MeV]; 
we take into account the uncertainty in the h b {lP) mass. A summary of the systematic 
uncertainties is given in Table IIV1 To obtain the total systematic uncertainty we add all 



TABLE IV: Systematic uncertainties in the rjb(lS) parameters from various sources. 





N[ Vb (lS)}, 10 3 


M[r] b (lS)}, MeV/c 2 


r[ Vb (lS)], MeV 


Range in M m i ss (7r + 7r ) fits 


+0.0 
-0.6 


+0.0 
-0.2 


+0.0 
-0.1 


Poly order in M m i ss (7r + 7r _ ) fits 


+0.0 
-0.6 


+0.1 
-0.1 


+0.0 
-0.4 


AM miss (7r + 7r~7) binning 


+0.2 
-0.1 


+0.3 
-0.8 


+1.0 
-0.8 


Range in AM m i ss (7r + 7r~7) fits 


+0.9 
-0.2 


+0.1 
-0.1 


+1.4 
-0.3 


AM m i ss (7r + 7r~7) bg parameterization 


+5.5 
-1.4 


+0.5 
-0.2 


+10.9 
-2.2 


Resolution function 


+0.2 
-0.1 


+0.5 
-0.5 


+1.4 
-1.4 


Selection requirements 




+0.4 

-1.9 


+3.0 
-2.0 


hb(lP) mass 




+1.1 
-1.1 




Total 


+5.6 
-1.7 


+1.4 
-2.4 


+11.5 
-3.4 



sources in quadrature. 

We study the shift of the rib(^S) parameters in case other line-shape parameteriza- 

: BW(m) EEo+ ( E 1 Eq )2 ; where 

BW(m) is the Breit-Wigner function, E [Eq] is the 7 energy in the h b (lP) rest frame 
[calculated for the pole mass of the r]b(lS)]. We also consider the CLEO parameteriza- 
tion [16]: BW{m) E 3 exp(— ™), where f3 is a fit parameter. Both the KEDR and CLEO 
Collaborations used these parameterizations for the J/tp — » ^7 transitions. We do not find 
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considerable shifts in the rjb(lS) yield, mass or width if these alternative parameterizations 
are used instead of the non-relativistic Breit-Wigner function. 

The significance of the rjb(lS) signal including systematic uncertainties is 13 a. 



MEASUREMENT OF B[h b (lP) r) b (lS)>y] 

We measure the rjb(lS) [/t&(lP)] yield in the events that fail the R2 < 0.3 or 7r° veto 
requirements [R 2 < 0.3 requirement]. In the fit to the AM miss (7r + 7r~7) [M miss (7r + 7r - )] dis- 
tribution for the rejected events we allow the rjb(lS) mass and width [the hb(lP) mass] to 
float within the uncertainties of the measured values given in Table IIHI (Table H]). We find 
Ni[r] b {lS)} = (5.5 ± 2.7 ± 2.1) x 10 3 and N^h^lP)] = (13.0 ± 2.8 ± 0.5) x 10 3 . 

From the Ni values and the yields from Tables |T] and IIHI we determine the total 
yields No[hb(lP)] and No[r]b(lS)] without requirements on R2 or a tt° veto. We obtain 
B = N [r]b(lS)]/No[hb(lP)}/ e, where e is the reconstruction efficiency of the radiative pho- 
ton, which is found from MC to be 74.1% with 2% systematic uncertainty (the MC statistical 
uncertainty is negligible). We find 

B[h b {lP) -> r]b(lS)j} = (49.8 ± 6.8^° 2 9 )%. 



CONCLUSIONS 



We report the first observation of the radiative transition hb{lP) — > 775(15)7, where the 
hb(lP) is produced in T(55) — > hb{lP)n + n^ dipion transitions. We report the single most 
precise measurement of the rjb(lS) mass, 

(9401.0±1.9l2l) MeV/c 2 , which corresponds to the 
hyperfine splitting AMm?\r)h(lS )] = (59.3±1.9j^; 4 ) MeV/c 2 . This value deviates significantly 
from the current world average 17] but decreases tension with theoretical expectations 
(see Fig. [7]). We report the first measurement of the rjb(lS) width (^A^l'^^f) Me V, which 
is in the middle of the range of predictions from potential models, 4 — 20 MeV (l8j. For 
the branchingfraction we find B[hb(lP) — > 776 (1 5)7] = (49.8 ± 6.8l5° 2 9 )% in agreement with 
expectations (^j. 

We also report updated results for the h b (\P) mass (9899.0 ± 0.4 ± 1.0)MeV/c 2 and 
hyper-fine splitting AMhf[^(1-P)] = (0.8 ± 1.1) MeV/c 2 . The latter is consistent with zero, 
as expected. 
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II • II 



I .... I 



BaBarY(3S) 



BaBarT(2S) 



CLEO 



Belle preliminary 



30 40 50 60 70 80 
AM HF (r| b ), MeV/c 2 

FIG. 7: Hyperfine splitting measured by BaBar in T(35) data [4], BaBar T(25) data [5], CLEO 
and present preliminary result of Belle. In addition, pQCD [?J (horizontally hatched) and Lattice 
QCD [8| (vertically hatched) predictions are shown. 
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